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Abstract. Electromagnetic localization and existence of gap solitons in nonlinear 
metamaterials, which exhibit a stop band in their linear spectral response, is 
theoretically investigated. For a self-focusing Kerr nonlinearity, the equation for the 
electric field envelope with carrier frequency in the stop band - where the magnetic 
permeability n(ui) is positive and the dielectric permittivity e(ui) is negative - is 
described by a nonlinear Klein-Gordon equation with a dispersive nonlinear term. A 
family of standing and moving localized waves for both electric and magnetic fields is 
found, and the role played by the nonlinear dispersive term on solitary wave stability 
is discussed. 
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1. Introduction 



The study of electromagnetic propagation in metamaterials, i.e. artificially constructed 
media showing negative effective dielectric permittivity e(uj) and magnetic permeability 
/Jj(u>), has received a tremendous and increasing interest in the past few years especially 
after the fabrication of microstructured materials showing negative refraction at 
microwaves [HO [3]. Electromagnetic wave propagation in left-handed (LF) materials 
in the linear regime has been extensively investigated since the pioneering work by 
Veselago [4] (see, e.g., [8] and references therein), and many unusual effects with 
a potential impact in different fields of science have been predicted, the most notably 
one being the realization of a "perfect" lens [9] which has raised a debate among the 
scientific community [10]. Only very few papers [TTJ [12] have been so far devoted to 
study nonlinear electromagnetic propagation in metamaterials despite they are expected 
to exhibit nontrivial nonlinear properties |13| . In particular, surface-polariton solitons 
have been studied in metamaterial interfaces assuming a simple Kerr-type nonlinearity 
for the dielectric response of the medium [12]. An interesting property of composite 
metamaterials, which is usually found in periodic media, is the existence of a frequency 
stop band in their linear spectral response corresponding to negative values of either 
e{uj) or h{uj). As compared to most common band gaps found in periodic media, where 
forbidden frequencies arise due to multiple Bragg scattering, in composite lossless LH 
materials forbidden frequencies exist because of the special dispersion curve for the 
effective permittivity and permeability, which have opposite sign in a spectral interval 
which separates allowed (propagative) regions of LH (i.e., e, fi < 0) and right-handed 
(RH, i.e. e, > 0) waves. Experimental evidence for such band gaps at microwaves 
was reported in Ref.pQ using a microstructured medium made of a periodic array of 
conducting split ring resonators and wires [14]. A rather universal feature of nonlinear 
periodic media exhibiting a frequency stop band in their linear dispersive response is the 
existence of self-transparent gap soliton envelopes supported by the nonlinearity [15J, 
which have been extensively investigated in the fields of nonlinear optics [15j [16j [171 [IB] 
and matter waves [T9l [20] and usually modelled by a generalized massive Thirring 
model. The existence of a frequency stop band in metamaterials connecting RH and 
LH propagative regions raises the question whether, in presence of nonlinear effects, gap 
solitons can be supported in microstructured media, despite Bragg scattering effects do 
not play any role and the medium behaves as a homogeneous one [H] . In this Letter it 
is shown indeed that envelope gap solitons can be supported in lossless metamaterials 
exhibiting a cubic dielectric nonlinearity when the carrier frequency of the fields is 
tuned inside the stop band separating LH and RH propagative spectral regions. In Sec. 
2 a nonlinear envelope equation for the electric field is derived starting from Maxwell's 
equations in the limit of a narrow spectral band gap. The resulting equation is a complex 
nonlinear Klein-Gordon equation (NLKGE) with a dispersive cubic term, which arises 
due the strong dispersion of the magnetic permeability for linear waves. Solitary wave 
solutions to this equation are studied analytically and numerically in Sec. 3. In particular, 
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a two-parameter family of steady or moving bright solitary waves for both electric and 
magnetic fields is found. Numerical simulations indicate also that these solitary waves 
are stable. 

2. Electromagnetic wave propagation in metamaterials with a cubic 
nonlinearity: derivation of the nonlinear envelope equation 

Let us consider propagation of electromagnetic plane waves in a dispersive medium, with 
frequency-dependent relative dielectric permittivity e{u) and magnetic permeability 
h(uj), which exhibits an instantaneous cubic (Kerr-type) nonlinearity in its dielectric 
response. Assuming quasi-monochromatic electric and magnetic fields with carrier 
frequency ujq polarized along the transverse x and y directions, 
E(z, t) = £(z, t) exp(iuot)u x + c.c. and 

H(z,t) = TC(z,t) exp(iuJot)u y + c.c, a nonlinear polarization term of the medium 
Pnl — eoX (E-E)E (x > for a self-focusing nonlinearity), from Maxwell's equation 
the following coupled-mode equations for the slowly-varying electric and magnetic field 
envelopes £ and TC can be derived in the rotating-wave approximation: 



where eo and /j,q are the vacuum dielectric permittivity and magnetic permeability, and 
the operators e(uj — id t ), [i{ujq — id t ) are defined, as usual, by the power expansion of 
e(uJo + Q) and /jl(ujo + Q) at around Q = with the substitution Q — > —id t (see, e.g., 
[2T]). A single equation for the electric field envelope E can be derived from Eqs.(l) and 
(2) and reads explicitly: 



where we have set k 2 {uj) = eo/iok-' 2 e( a; )M u; )- For a composite metamaterial made by 
an array of conducting nonmagnetic split-ring resonators and continuous wires, we may 
assume the following general form for the permeability and permittivity (see, e.g. [1] 
and references therein): 



where oj p , u m , F, 7 e and 7^ can be tuned by changing the geometrical and physical 
parameters of the microstructures forming the composite medium [H [3] . For an ideal 
lossless medium, one has % = 7^ = 0. Referring to this ideal case and assuming a plasma 
frequency lu p larger than ZJ = u m /(l — F) 1 / 2 , linear waves of Eq.(3) are evanescent 
{k 2 {uj) < 0) in the spectral interval uJ < uj < u p (band gap) , where ji{oj) > but 
e(u) < 0, whereas they are propagative {k 2 {uj) > 0) for uj > uj v and uj < ZJ (down to the 
resonance uj m ) , with a LH [RH] behavior of the material on the left [right] side of the band 
gap. A typical behavior of e(uj) and fi(uj), for parameter values compatible with recent 
experimentally fabricated structures p], is shown in Fig.l. Let us now assume a reference 




(1) 
(2) 



d 2 £ + k 2 (uj - idt)£ = -3^ eo/ioX (3) (^o - id t )n(coo - id t )\£\ 2 £, 



(3) 
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carrier frequency u> for the fields at the center of the band gap, i.e. uj = (u p +TJ)/2, and 
assume a sufficiently narrow band gap {oj p — To <C lo ) such that the dispersion curves e{uj) 
and ii(uj) can be expanded up to leading order in lo— ojq at around the bandgap region, i.e. 
e(io) ~ e' (LO — LO — A) and fx(to) ~ — loq + A), where A = (lo p — To)/ 2 is the half width 
of the stop band and e' = (de/du}) Uo , /i = (d/i/dco)^ are real- valued for the lossless 
medium and always positive for causality [22]. Such approximation corresponds to 
assume a parabolic behavior of k 2 {uj) versus lo, which is reasonable as shown in Fig. 1(b). 
With these assumptions, for slowly-varying envelopes ( (l/u>o)dt ~ (A/uq) <C 1) and 
introducing the dimensionless variables t' = At, z' = Az/v g , ip = [3x^ /(g'qA^^E, and 
(p = [3x (3) yu / '{e e' A)] l ^H, where v g = {e fi io%e ii' )- 1/2 , Eqs.(l) and (2) read: 



d z xp 




•1 + idf) <p 



1 +id t r) ip — i* 



\ Mo 



(5) 



(6) 



Finally, from Eqs.(5) and (6) the following nonlinear equation for the normalized electric 
field envelope ip can be obtained: 



1 



ip = id t < {\ip^ip 



(7) 



which is the basic equation describing nonlinear wave propagation in Kerr-type nonlinear 
metamaterials at frequencies close to (or inside) the stop band separating the LH and 
RH propagative spectral regions of the medium. 



3. Gap solitons 

In this section we focus our analysis on the existence and stability of stationary and 
moving solitary waves to the nonlinear envelope equation derived in the previous section, 
Eq.(7). From a mathematical viewpoint, Eq.(7) is a complex NLKGE which differs from 
the most usual ones encountered in other physical fields (see, e.g. [23J, [2H [251 126] ) 
due to the presence of a nonlinear dispersive term, represented by the right hand 
side term in Eq.(7). The physical reason for the appearance of a nonlinear dispersive 
term in the envelope equation despite the instantaneous (i.e. non-dispersive) nature 
of the Kerr nonlinearity is due to the fact that close to the bandgap (u ~ u>o) the 
magnetic permeability fJ>(u) is small and its frequency dependence can not be neglected, 
leading to the appearance of the nonlinear dispersive term in Eq.(7) [27]. Using a 
common terminology adopted for solitary waves in nonlinear periodic media, we will 
call stationary or moving localized waves to Eq.(7) gap solitions. However, the reader 
should be aware that gap solitons in metamaterials, as described by Eq.(7), have 
a different physical origin and should not be confused with usual Bragg solitons in 
nonlinear periodic media, which are described by a generalized massive Thirring model 
for counterpropagating fields. In fact, as we already pointed out [H], despite artificial 
metamaterials show a microscopic periodic structure, they behave as homogeneous 



S. Longhi, Gap solitons in metamaterials 



5 



media and Bragg scattering does not occur. In this sense, gap solitons to Eq.(7) should 
be regarded as a kind of self-transparency solitary waves rather than true Bragg solitons, 
thought no coherent effects of the medium are involved. 

In order to find an analytical form of gap solitons, let us first observe that, if the 
dispersive nonlinear term were negligible, Eq.(7) would reduce to the well-known if) A 
NLKGE; existence and stability of bright solitary waves for such an equation were 
previously studied in Ref.[23| l2lj . Steady solitary waves read explicitly 
ip(z',t') = [2(1 — fi 2 )] 1 / 2 sech[(l — fi 2 ) 1 / 2 ^'] exp(zfW'), where f2 is a free-family parameter 
(\Q\ < 1 for existence), which measures the frequency offset of the soliton from the 
band gap center. By exploiting the Lorentz invariance of the ip A NLKGE, a family of 
moving gap solitons, with an arbitrary velocity v satisfying the condition \v\ < 1 (slow 
gap solitons), can be then generated. A linear stability analysis of such solitary waves 
shows that they are linearly stable provided that l/y/2 < |Q| < 1 [24J. If -as it is our 
case - the nonlinear dispersive term in Eq.(7) can not be neglected, an explicit analytical 
form of solitary waves can not be obtained in general, and Lorentz invariance is broken. 
We can nevertheless look for moving bright solitary waves of Eq.(7) in the form [28] 
if}(z', t') = F(z' — vt') exp[i(Qt' — Qz')] with F(£) — > as £ — > oo, where v is the solitary 
wave speed and Q, Q are real-valued parameters; with the further constraint Q = vfl, 
the envelope F can be found as a homoclinic loop emanating from the unstable solution 
F = of the equation: 

7^ - (1 - 1&)F + ( 1 + Q + iv-^) \F\ 2 F = 0, (8) 



where we have set 7 = 1 — v 2 . An inspection of the asymptotic (linear) form of Eq.(8) as 
£ — > 00 shows that solitary waves exist provided that the condition (1 — 7^ 2 )/7 > is 
satisfied, which implies 7 > 0, i.e. \v\ < 1 (slow gap solitons) and \ fl\ < l/y/j- For fixed 
values of the two family parameters Q and v , the homoclinic trajectory of Eq.(8) can be 
numerically computed by standard techniques; the corresponding magnetic field tp can 
be then computed by numerical integration of Eq.(5). Explicit analytical expressions of 
companion electric and magnetic gap solitons have been found solely in the steady case 
(v — 0), which read: 

tP(z', t') = [2(1 - fi)] 1 / 2 sech[(l - n 2 ) l ' 2 z'} exp(^t') (9) 



y/2(l-Q) sinh[(l-0 



2\l/2 



z' 



= - H i Jrm coA-ffWi exp(iSi(,) - (10) 

Note that the steady solitary waves (9) and (10) describe a stationary nonlinear 
localized mode for the electromagnetic field provided that its carrier frequency falls 
inside the bandgap region (\Q\ < 1). Examples of electric and magnetic field profiles, 
for a stationary [v = 0; see Eqs.(9) and (10)] and moving (v 7^ 0) gap solitons, are 
shown in Fig. 2. Note that, in terms of real physical variables, spatial and temporal 
length scales in the figures, for parameter values as in Fig.l, are 5.77 cm and 0.32 ns, 
respectively. The soliton velocity v is measured in terms of the characteristic velocity 
v g = (u;Qeo/ioeoA i o)~ 1//2 ' whose value is v g ~ 0.6cq, where Co = l/(eo/io)^ 2 is the speed of 
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light in vacuum. 

An important issue is the stability of the solitary waves. A rigorous stability analysis is 
challenging and can not be framed in the analysis of Klein-Gordon solitons developed 
in Ref.[23]. However, from numerical integration of Eqs.(5) and (6) we could ascertain 
stable propagation in a wide range of the existence domain. In particular, we found that 
the dispersive nonlinear term entering in Eq.(7) plays a stabilizing effect on the solitary 
wave dynamics. We numerically integrated Eqs.(5) and (6) using a pseudospectral 
split-step technique with typical 512 discretization points in a traveling reference frame 
£ = z' — vt' and rj = t', where the unperturbed solitary wave is at rest; as an initial 
condition we assumed a solitary wave perturbed with a small random noise to seed 
possible instabilities. As an example, stable evolution of a steady solitary wave, with 
carrier frequency at the band gap center, and of a moving solitary wave are shown in 
Fig. 3. We note that, if the nonlinear dispersive term in Eq.(7) were neglected, the 
steady solitary wave of Fig.3(a) would be unstable according to the analysis of Ref.[24j. 
We checked indeed the stabilizing effect played by the nonlinear dispersive term by a 
numerical analysis of Eq.(7) with and without the right hand side term, assuming as an 
initial condition their respective non-moving solitary waves with an added small random 
noise. We assessed the reliability of our numerical method by exactly reproducing the 
stability domain l/y/2 < \Q\ < 1 for steady solitary waves of the V' 4 NLKGE according 
to the linear stability analysis [21]. Figure 4 shows the numerical results of evolution 
of the steady solitary waves of Eq.(7) at the band gap center in absence [Fig. 4(a)] and 
in presence [Fig.4(b)] of the nonlinear dispersive term. Note that, as in the former 
case the solitary wave is unstable according to the linear stability analysis [24j, in 
the latter case no unstable growing modes were observed in the numerical simulation, 
indicating that the nonlinear dispersive term plays a stabilizing effect. The emergence of 
instabilities for gap solitons and their physical explanation are usually nontrivial issues 
and represent a challenging task |29j; therefore a detailed physical explanation of the 
instability suppression observed in our numerical simulations goes beyond the aim of the 
present work. We just mention that for the NLKGE without the nonlinear dispersive 
term the instability of the gap soliton close to the band gap center < l/y/2), 
as numerically reproduced in Fig.4(a), arises because of the emergence of an unstable 
internal mode of the soliton [21]. Our numerical results suggest that the presence of 
the nonlinear dispersive term in the NLKGE makes such internal mode damped in the 
entire bandgap region of linear waves < 1). 

4. Conclusions 

In conclusion, we have predicted the existence of stationary and moving gap solitary 
waves for electric and magnetic fields in metamaterials with a Kerr nonlinearity. These 
waves, which are supported by the nonlinear dielectric response of the medium, exist 
in the band gap spectral region of the medium, corresponding to e < and fi > 0, 
which separates the spectral regions of allowed propagation where the medium behaves 
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as a LH or as a RH material. The gap solitary waves studied in this work - which 
can be described by a complex NLKGE with a dispersive cubic term - represent a new 
class of gap solitons in artificially-constructed microstructured periodic media which do 
not involve Bragg scattering and provide a noteworthy and physically relevant example 
of self-transparent electromagnetic wave propagation in the recently developed class of 
composite metamaterials. 
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Figure 1. Typical dispersion curves of a composite metamaterials with a stop band 
in the microwaves, (a) Behavior of relative dielectric permittivity e(u>) and magnetic 
permeability /u(w) versus frequency in the ideal lossless medium for parameter values 
lo p = 2ir x 10 GHz, uj m = 2tt x 6 GHz and F = 0.56 (corresponding to Zu ~ 2ir x 9 
GHz and a central bandgap frequency ujq ~ 2ir x 9.5 GHz). The dashed curves 
(almost overlapped) in the figure show, for comparison, the behavior of imaginary 
parts of e and fi in a low-loss medium with damping terms 7 e /w p = j^/wp — 0.005. 
(b) Corresponding behavior of k 2 (ui) versus frequency for linear waves in the lossless 
medium. 
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Figure 2. Gap solitary waves for normalized electric field tp (upper figures) and 
magnetic field tp (lower figures) for a few values of family parameters v and Q. 
Continuous curves: modulus square of the fields ; dotted curves: real part of the 
fields; dashed curves: imaginary part of the fields. In (a) steady soliton with v = 
and = (see Eqs.(9) and (10) given in the text); in (b) moving soliton with v = 0.8 
and f2 = 0; in (c) moving soliton with v = 0.8 and SI = 1.2 (the real and imaginary 
parts are taken at t' = in this case). In the plots of lower figures (magnetic field), 
we assumed e /fi Q = 1. 
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Figure 3. Numerically computed evolution of gap solitons of Eqs.(5) and (6) in the 
moving reference frame £ = z' — vt', r\ = t' for (a) v = = and (b) v = 0.4, = 0.5. 
Discretization points: 512. Time step: dr\ = 0.002. Initial condition is the solitary 
wave with an added small random noise (1% amplitude). 
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Figure 4. (a) Temporal evolution of the stationary solitary wave of the ip 4 complex 
NLKGE [Eq.(7) with right hand side equal to zero] for O = (time step dt' — 0.001, 
512 discretization points), (b) Same as in (a), but for Eq.(7). 



